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It is well known that the particles in a beam of Boson obey- 
ing Bose-Einstein statistics tend to cluster (bunching effect), 
while the particles in a degenerate beam of Fermion obeying 
Fermi-Dirac statistics expel each other (anti- bunching effect). 
Here we investigate, for the first time, the statistical correla- 
tion effect for the composite Boson, which is formed from a 
spin singlet entangled electron pair. By using nonequilibrium 
Green's function technique, we obtain a positive cross corre- 
lation for this kind of the composite Boson when the external 
voltage is smaller than the gap energy, which demonstrates 
that a spin singlet entangled electron pair looks like a com- 
posite Boson. In the larger voltage limit, the cross correlation 
becomes negative due to the contribution of the quasiparti- 
cles. At large voltages, the oscillation between Fermionic and 
Bosonic behavior of cross correlation is also observed in the 
strong coupling regime as one changes the position of the reso- 
nant levels. Our result can be easily tested in a three-terminal 
normal-superconductor-superconductor (N-S-S) hybrid meso- 
scopic system. 
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I. INTRODUCTION 

There are two kinds of quantum statistics in nature. 
All particles have either half-integral or integral spin (in 
units _pf the Plank constant Ti) and they obey Fermi- 
DiracEl or Boson-Einsteino statistics respectively. It is 
also noteafl that there is an effective attraction be- 
tween the Bosons and an effective repulsion between the 
Fermions. This is the well known statistical correlation 
effectso, which are purely quantum effect. The experi- 
ments examining the quantum statistical properties date 
back to. the pioneer work, by Hanbury Brown and Twiss 
(HBT)lj. They used photon intensity interferometry to 
probe the intensity correlation information between two 
partial beams, which was generated by a beam splitter. 
Due to the Bosonic property of photon, the positive in- 
tensity correlation was observed, indicating an enhanced 
probability for the simultaneous detection of two pho- 
tons, one in each partial beam. This means that photons 
tend to bunch in cluster. Several theoretical works have 
suggested the different analogies of this experiment with 
electrons in mesoscopic systems. The Fstmionic analog 
of HBT experiments, one by Henny et al.B and the other 
by Oliver et al.Q, showed the expected negative inten- 
sity correlation and observed the anti-bunching effect. 
In this paper we will investigate the HBT experiment for 
the composite Boson. This composite Boson is formed 
by a spin singlet entangled electron pair, which will be 



discussed below. Due to the zero total spin for the en- 
tangled electron pair, we expect that these composite 
Bosons tend to bunch in cluster. Motivated by the ap- 
plication in quantum communication and computation, 
Burkard et alEl have studied entangled electrons in an in- 
teracting many-body environment. They found that the 
Fano factor for singlets is twice as large as for indepen- 
dent classical particles and is reduced to zero for triplets. 
Torries and Martina investigated a three-terminal N-N-S 
mesoscopic system, both positive and negative correla- 
tions were found in the Aadreev regime. Very recently, 
Samueisson and ButtikertH studied the same structure 
and found the positive correlation for a wide range of 
junction parameters which survives even in the absence 
of the proximity effect. The statistics of charge trans- 
port of a three terminal N-N-S beam splitter has also 
been investigatecLJ and positive cross correlation is found 
between the currents in two normal leads for a wide pa- 
rameter range. Instead of the structures of the references 
here we consider a three-terminal mesoscopic N-S-S 
hybrid system. This structure is a direct photon analogy 
of the HBT interferometer which has a normal lead and 
two superconducting leads. A quantum dot, connected 
by these three terminals, acts as a splitter. Suppose 
that the chemical potentials fi s for both superconduct- 
ing reservoirs are set to zero, and the chemical potential 
for the normal is above zero, i.e., eV > 0, which guaran- 
tees the electron current passing from the normal lead to 
both superconducting leads. We further assume the tem- 
perature is very low. If the external voltage eV is smaller 
than the gap energy A of the superconducting leads, the 
single quasiparticle current is forbidden. In this case, 
we only have two-electron current due to the presence of 
Andreev reflection process, i.e., incoming electrons being 
Andreev reflected into outgoing holes with the transfer 
of a Cooper pair into the superconductor. This means 
that an electron (with energy e above the Fermi level 
and spin a) in the normal lead, has to combine another 
electron (with energy — e, below the Fermi level and spin 
—a) to pass through_|thc NS interface. Does this en- 
tangled electron paiiO look like a composite Boson? or 
rather, can we obtain a positive cross correlation function 
(<AI a AIf} > with between two superconducting 

leads? The purpose of this paper is to answer this ques- 
tion. We note that due to the current conservation the 
cross correlation function of a two-lead system must be 
negative regardless of normal or superconducting leads. 
Instead of considering the fluctuation in a single electron 
beam through the two-lead system, the HBT experiment 
considered here focuses on the cross correlation of two 
beams from the beam splitter. Hence we expect positive 
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cross correlation at small voltages which is indeed what 
we found in this work. When eV > A, the quasiparticles 
will also participate the transport. Due to the Fermionic 
nature of quasiparticles, it will partially cancel the pos- 
itive contribution of the entangled electron pair to the 
cross correlation. The competition of these two contribu- 
tions from the entangled electron pair and quasiparticles 
can lead to either positive or negative cross correlation 
depending on which contribution dominates. 



II. THEORETICAL FORMULATION 

We begin with the following model Hamiltonian 

H = X/ € P C l,p<7 G l,p<7 + X/E 6 k C nM Cn M 
p kn a 

+AC+ feT C+_ fci + AC n ^ kl C nM ] + eo^X 

+E[ T ^ c iV^ + c - c -i + I2i T ^c+ ka d„ + c.c] (i) 



where the first term denotes the Hamiltonian of the nor- 
mal lead. The second term (n — 2, 3) describes the 
Hamiltonian of two BCS superconducting leads. Here 
C\ u a is the creation operator of electrons in the nor- 
mal lead and G^ fe is the corresponding creation opera- 
tor in the superconducting lead. The third term is the 
Hamiltonian for a quantum dot, which is used to mimic 
a tunable beam splitter. Here we have applied a gate 
voltage which can control the level of the dot so that 
e = e^+evg. Without loss of generality, we set 6q = 0. 
The other terms in Eq.([l]) are Hamiltonians describing 
the couplings between the quantum dot and leads. To 
simplify the discussion, we have assumed that two su- 
perconducting leads have the same gap energy A. We 
have also neglected the supercurrent between two super- 
conducting leads£3 and assumed that the hopping matrix 
elements are independent of the spin index. 

In the following, we will calculate the cross correlation 
between two partial beams through two superconducting 
leads. The current operator for the superconducting lead 
2 or 3 is 



I a = I a t(t) + I a i(t) 



with 



The cross correlation between two superconducting leads 
is defined as 

p 23 =< Ai 2 (h)Ai 3 (t 2 ) >=< [i 2 {h) - h][h{t 2 ) - h] > 

with I a =< I a >. Here < ... > denotes both the statisti- 
cal average and quantum average on the nonequilibrium 
state. Using the expression of the current operator, the 
cross correlation between two superconducting leads is 

P 23 = -4e 2 ^[T 2 , fc T 3 ^G< m (2,l)G> fc , T (l,2) 

k,k' 



-T , 2 i feT 3 * fe ,G^ TfeT (2, l)G^ d j(l, 2) 
_ ^2*fe^3,fe'G^ TdT (2, l)G^ Tfe , T (l, 2)] 



(2) 



where we have used abbreviation G(ti,t 2 ) — G(l, 2) and 
we have used k and k' to label, respectively, the second 
and third superconducting lead. The Green's functions 
Qr : a : < : > ixi 9_x 2 Nambu representation take the follow- 
ing formsEJcj 



<{X a1 (t 1 ),Y+(t 2 )}> <{X a i(t 1 ),Yp l (t 2 )}> 
< {*+ (h ) , Y+ (t 2 )} >< {X+ (tO, Y PI (t 2 )} > 



r < ( . , x _ ■ (<Y+(t 2 )X a ^(t 1 )> <Y f}i {t 2 )X a] {t l )> 
ap{ l! 2j l {<Y+(t 2 )X+( tl )> <Y Pl (t 2 )X+( tl )> 



<x a1 (h)Y+{t 2 )> <x QT (t 1 )r /3i (t 2 )> 

<Kl (*i)^3 + T <Ki (h)Y fu (i a )> 



where X and Y stand for the annihilation operators, such 
as Gi, p , C nt k, and d. These Green's functions satisfy the 
general motion, G> = G<+G' G a . Using the Keldysh 
equatiorO 

G<'> = (1 + G r S r )G^' > (l + S Q G a ) + G r S<G a 
we have the following relations 



^W^ii^) - T tk J dt [G r d ^ d i i (ti 1 t)g k ^ k ^(t,t 2 ) 

+Gdidi (*i i *)5fctfeT (*> *2) + G r d ^ di {t\ , t)g k ^ (t,t 2 ) 
+G d i dl (h,t)g a kW (t,t 2 )} 



(3) 



where a = 2, 3. Due to the electron-hole symmetry of 
the system, we have I a ]{t) = Iai(t). Hence the current 
operator can be rewritten as 

I a (t) = 2ie}}T a kCg_ kr di - c.c] 



G k j d ^(ti,t 2 ) — T 2 ,k J dt [5 feT : feT (ti, t)G d]d1i (t,t 2 ) 
+9 r km (h,t)G<i>(t,t 2 )] 



(4) 
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+G<i> l (t 1 ,t)g a k , lk ,^t,t 2 )] (5) 
where G k * da is given by 

Cfc T(icr (ii,t2) = T^h J dt[gl 1k ^(ti,t)G r d1drJ (t, t 2 ) 

+g r km (h,t)G dlda (t,t 2 )] (6) 

Substituting the above relations to Eq.(l) and taking 
the Fourier transform, we obtain 

P 23 = -4e 2 r 2 r 3 J g{(G I 'g < +G < g a ) TT (G r g> 

+G>g a ) TT + (g<G a + g r G<) TT (g>G a + g r G>) TT 
-G> {(g r G'-g<) TT + [(g<G a + g r G<)g a ] TT } 



-G T < T {(g r G'-g>) Tt 



[(g>G Q +g''G>)g Q ] TT }} (7) 



where r a = 27r PNa\T a k\ 2 with a = 2,3 are the 
linewidth functions. Here PN2,3 are the normal density 
of states of the superconducting leads 2 and 3. We have 
used the wide-band limittU and thus the linewidth func- 
tion is independent of the energy. G%£} < ' > = G r d ^ d ^i > are 
the full Green's functions for the quantum dot in pres- 
ence of the leads, while <j£'"/ are the exact Green's 
functions for BCS superconductor in the absence of the 
coupling between the leads and quantum dot. Eq.(Q) is 
the central result of this paper. It describes the cross cor- 
relation for a three-terminal hybrid N-S-S system and is 
valid at any temperature and finite voltage, i.e., valid for 
both eV > A and eV < A. In order to calculate this cor- 
relation, one must know all the Green's functions. The 
exact Green's functions for the isolated supercon- 

ducting leads araiSEJ 



f{E) 



- A 2 



E A 
A E 



g<(E) = if(E)9(\E\-A) 



VE 2 - A 2 



E A 
A E 



where f(E) = l/[exp(/3(_E — Ep) + 1] is the well known 
Fermi distribution function, 9(x) is the step function and 
((E) = 1 when E > -A, otherwise ((E) = -1. We will 
choose the Fermi energy of the normal lead in line with 
the chemical potential fi s of superconducting condensate 
which is set to zero, i.e., Ep = /j, s = 0. The retarded 
Green's function for the quantum dot can be calculated 
using the Dyson equation 



G r (E) 



1 



G r (E) 



1 



E-e 
E + e 



and 



1 



((E) 



E A 



The lesser Green's function can be obtained from the 
Keldysh equation G < = G r S < G a . Here the lesser self- 
energy is given by 



£<(£) = iTt 



ff(E + eV) 







f(E-eV) 



if(E)9(\E\-A)^0^((E) (I* 



Let's first consider the case that external voltage is 
smaller than the gap energy and consider zero tempera- 
ture behavior so that there is no quasiparticles partici- 
pating the transport. In this case, only two-electron cur- 
rent exists, i.e., the current from incoming electron and 
Andreev reflected hole, we have g r = g a and g <!> = 0. 
using the fact that 



G ,T,G r ( f+ ,° ! G" 



and 



G> = T\G r 



r (f+ - 1 



/_ - 1 
Eq.(^) can be further simplified as 



_ e 2 r 2 r 2 r 3 A 2 r dE 
p 23- A2 _ E2 j 

\/~ir r~ia s~ir r~ia |2 



(8) 



(9) 



4e 2 r 2 r 3 



(r 2 + r 3 

xT A (E)(l-T A (E)) 



/g,- ( w + , 



(10) 



where Ta(E) = TlG^G^ is the Andreev reflection co- 
efficient and f±(E) = f(E ± eV). Just as we expected, 
Eq.(|To|) is a positive quantity. To get more physical in- 
sight, we will assume that eV are small enough and we 
will keep only the first order in V in Eq.(fio|). We have 



p2Z — 



r 2 r 2 r 3 e 3 y 



, e 2 [ r 2 + (r 2 + r 3 ) 2 ] [r 2 - (r 2 + r 3 ) 2 ] 2 

2 ! 4 



>{4 + 



} (ii) 



For eV > A, we have to calculated P23 numerically which 
is presented in the next section. 



with 
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III. RESULT AND DISCUSSION 



We first use Eq.(lO) to calculate the cross correlation 
at finite voltage while keeping eV < A. In the follow- 
ing, we will use A as the unit of energy and study the 
symmetric case where r 2 = IV In Fig.l we show the 
cross correlation versus the gate voltage at fixed external 
bias eV = 0.6. Four different sets of coupling constants 
r are chosen: (1). Ti = 0.8 and r 2 = 0.8. In this case, 
it represents the strong coupling between leads and the 
quantum dot. The cross correlation (solid line) displays 
two broad peaks located symmetrically at ev g ~ ±0.6. 
(2). Ti = 0.8 and r 2 = 0.1. In this case, the normal 
lead couples strongly with the quantum dot while the 
superconducting leads couple weakly. The cross correla- 
tion (dotted line) exhibits a single peak at v g — 0. (3). 
Ti = 0.1 and r 2 — 0.8. This is the reverse of case (2). 
We see that the cross correlation (dot-dashed line) shows 
a flat region near v g = 0. (4). Ti =0.1 and T 2 = 0.1. 
This is the weak coupling case. The cross correlation 
(dashed line) has two sharp peaks close to v g = and 
decays quickly away from it. To understand these fea- 
tures, we notice that two terms (Fx = J dE Ta{E) and 
F 2 = JdE T\(E)) in Eq.@ tend to cancel each other. 
In the strong coupling case, the contribution from both 
terms are of the same order of magnitude. Both show 
broad peak near v g — with the second term decreasing 
faster away from v g = (see Fig. 2a). As a result, we ob- 
tain the double-peak structure as shown in Fig.l. When 
Ti = 0.8 and T2 — 0.1, the contribution from F% is much 
smaller than that of F\ and hence just one peak shows 
up (Fig.2b). For Ti = 0.1 and T 2 = 0.8, the addition 
of Fi and F 2 gives a long plateau between ev g = —0.5 
and ev g — 0.5 (Fig. 2c). In the weak coupling regime, Fi 
and F 2 give comparable contributions with a single peak. 
Similar to the strong coupling case, the integral of F2 
decreases faster than that of F\ resulting again a double- 
peak structure (Fig. 2d). As one decreases the external 
bias to eV — 0.2, the cross correlation in the strong cou- 
pling case still shows double peak structure with smaller 
amplitude (Fig. 3). We also find that the peak position is 
shifted towards origin and peak to valley ratio becomes 
much smaller. In the weak coupling regime, the cross 
correlation is roughly unchanged. In the case of Ti = 0.8 
and r 2 = 0.1 (or vice versa), the cross correlation de- 
creases. Now we examine the cross correlation versus 
external bias at fixed energy levels. Fig. 4a displays the 
cross correlation P 23 versus external voltage when e = 0. 
We see that, except for Ti = 0.8 and T 2 = 0.1 that P23 
increases monotonically, P23 develops a plateau region 
for the other three sets of coupling parameters. These 
plateau regions are due to the resonant tunneling which 
can be seen from Fig. 4b where the differential cross corre- 
lation dP2s/dV versus external voltage isjiepicted. Here 
we see typical behavior of the shot noisecd: a minimum 
separated by two peaks. The minimum is due to the reso- 
nant Andreev reflection since dP23/dV ~ Ta(1 — Ta). As 



one increases the energy level (eo = 0.3), the dip between 
two peaks can no longer reach zero indicating that the 
maximum Andreev reflection coefficient Ta is much less 
than one. We also note that for Ti = 0.1 and T 2 = 0.8, 
only one peak is left and resonant feature disappeared. 

To study the effect of quasiparticle when eV > A, we 
calculate the cross correlation using Eq.(j^). Fig. 6 shows 
the cross correlation versus external voltage at eo = 0.0. 
We see that once the voltage is larger than the gap en- 
ergy A, the cross correlation decreases quickly indicat- 
ing Fermionic contributions. For the strong coupling 
case, P 2 3 becomes negative in the large V limit. This 
can be understood as follows. When eV > A, electrons 
with energy less than eV will all participate in transport. 
In particular, for incoming electrons with energy inside 
the superconducting gap, only two-electron current is al- 
lowed and hence the contribution to the cross correlation 
should be positive as we just discussed above. However, 
when the energy of incoming electron is out side of the gap 
the current comes from of four processescjtLa: (1). An- 
dreev reflection; (2) . the conventional electron tunneling 
through the system; (3). " Branch crossing" processes: an 
electron incident from the normal lead converting into a 
hole like in the superconducting leads; (4). An electron 
(or a hole) incident from the normal lead tunnels into 
the superconducting lead, picks up a quasiparticle (or a 
quasihole) in the superconducting lead and creates (or 
annihilates) a Cooper pair. In these processes, the lat- 
ter three give negative contributions to the cross corre- 
lation. Competition between Andreev reflection process 
and the rest of three processes give rise either positive 
or negative cross correlation depending on which process 
dominates (see Fig. 6). Typically, near the resonance the 
Brei tp^ijj mer form for the Andreev reflection coefficient 



r 2 r 2 
1 i 1 2 



4(£ 2 - el + rjr/4) 2 + r 2 r 2 + e 2 (r + sry 



(12) 



and transmission coefficient for normal tunneling process 



T = 



rir 2 



(£-e ) 2 +r 2 /4 



(13) 



where T = Ti + T 2 and 6T = T x - T 2 . We see that 
the Andreev reflection is suppressed when off resonance. 
Furthermore, at large external bias if the resonant energy 
is outside of the gap, the Andreev reflection is drastically 
suppressed and normal tunneling is allowed at certain en- 
ergy. Therefore, we expect negative cross correlation in 
this case. In Fig. 7, we depict P23 versus V at e = 2.0. 
Since the resonant level is outside of the gap, the plateau 
region for P23 when eV is inside the gap disappeared. We 
see that except for the case of Ti = 0.8 and T 2 = 0.1, P23 
becomes negative at large voltages. Our numerical result 
shows that at even larger eo, the transport of quasiparti- 
cle dominates and all P23 are negative at large external 
voltage. Finally, we plot in Fig. 8 the P23 versus v g at 
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eV = 4. We see that at large voltages, all the cross cor- 
relation functions become negative. For the strong cou- 
pling case, we observe oscillations of P23 between Bosonic 
and Fermionic behaviors due to the competition between 
entangled electron pair and the quasiparticles. This can 
be easily checked experimentally by changing the gate 
voltage. 

In summary, we have proposed an entangled electron 
pair HBT experiment by using the three-terminal N-S-S 
hybrid mesoscopic system. When the external voltage 
is less than the gap energy, only two-electron current 
is present. The cross correlation is found to be pos- 
itive, which demonstrates that this entangled electron 
pair looks like the composite Boson and tends to bunch 
in cluster. However, when the external voltage is larger 
than the gap energy the quasiparticle will participate the 
transport which gives the Fermionic contribution to the 
cross correlation. As the result of competition between 
Andreev reflection process and the other tunneling pro- 
cess involving quasiparticles, the cross correlation can be 
either positive or negative depending on which one domi- 
nates. For the strong coupling case and at large external 
voltage, the cross correlation function changes sign as one 
varys the gate voltage which controls the position of the 
resonant level. 
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FIG. 2. The contribution of Fi (dotted line) and F 2 
(dot-dashed line) to the cross correlation P23 (solid line) ver- 
sus gate voltage at eV = 0.6. (a). Ti = 0.8 and T 2 = 0.8; 
(b). Ti = 0.8 and T 2 = 0.1; (c). Ti = 0.1 and T 2 = 0.8; (d). 
ri = 0.1 and T 2 = 0.1. 

FIG. 3. The cross correlation P23 versus gate voltage at 
eV = 0.2. The coupling parameters and corresponding sym- 
bols are the same as Fig.l. 
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FIG. 4. (a). The cross correlation versus external voltage 
at eo = 0.0. (b). The differential cross correlation versus 
external voltage at eo = 0.0. The coupling parameters and 
corresponding symbols are the same as Fig.l. 



FIG. 5. The differential cross correlation versus external 
voltage at eo = 0.3. The coupling parameters and correspond- 
ing symbols are the same as Fig.l. 



FIG. 6. The cross correlation versus external voltage at 
eo = 0.0. The coupling parameters and corresponding sym- 
bols are the same as Fig.l. 



FIG. 7. The cross correlation versus external voltage at 
eo = 2.0. The coupling parameters and corresponding sym- 
bols are the same as Fig.l. For illustration purpose, we have 
multiplied the cross correlation by a factor of 10 for the dotted 
line, 5 for the dot-dashed line, and 50 for the dashed line. 
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FIG. 8. The cross correlation versus gate voltage at 
eV — 4.0. The coupling parameters and corresponding sym- 
bols are the same as Fig.l. 
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